Ferromagnetic Polarons, Phase Separation, Stripes and Polaron Lattice: 
The two- and three-dimensional ferromagnetic Kondo Model 
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We investigate the two- and three-dimensional ferromagnetic Kondo lattice model by unbiased 
Monte Carlo simulations. A phase diagram for the two-dimensional model is presented, in which 
the stability of magnetic order and ferromagnetic polarons is examined with respect to the antifer- 
romagnetic superexchange J' and temperature. The Monte Carlo simulations reveal that J' > 0.02 
strengthens individual polarons while small J' < 0.02 favors larger clusters and phase separation 
except for small doping. Lowering the temperature stabilizes ferromagnetic polarons for realistic 
J' > 0.01, while phase separation is only favored for very small J' < 0.01. Our Monte Carlo 
simulations show that low temperatures can lead to diagonal or vertical stripes depending on J'. 

Simulations for three-dimensional systems yield ferromagnetic polarons, which form a 'polaron 
lattice' at higher doping levels 0.2 < x < 0.23, when independent polarons do no longer fit into the 
system. No tendency to phase separation is observed in three dimensions. 



I. INTRODUCTION 

Doped and undoped manganese oxides, such as 
Ri-^A^MnOa, where R denotes a rare earth and A an 
alkaline earth, have a very complex phase diagram 1,2 
depending on temperature, doping, ionic radius of the 
involved elements or magnetization. They have been 
thoroughly investigated because of the CMR effect and 
have been found to have a rich phase diagram includ- 
ing ferromagnetic (FM), antiferromagnetic (AFM), insu- 
lating, metallic and charge ordered states. Quasi two- 
dimensional (2D) systems with well separated MnOV 
(bi)layers exist. 

A comprehensive understanding of their properties cer- 
tainly requires treatment of the phononic degrees of free- 
dom and of the orbital degeneracy including Coulomb 
repulsion, which however would pose overwhelming diffi- 
culties, especially in more than one dimension. 

The much simpler FM Rondo lattice model 3 captures 
some of the properties of these materials, namely the hop- 
ping of the itinerant manganese e g electrons via the inter- 
mediate oxygen ions (double exchange) and their strong 
Hund's rule coupling to the S = 3/2 corespin formed by 
the localized manganese ti g electrons. An additional an- 
tiferromagnetic coupling between the corespins is often 
included to account for superexchange of the t?, g elec- 
trons. Even for this model, full quantum mechanical cal- 
culations are difficult in more than one dimension, and it 
has therefore been proposed to replace the ti g corespin by 
a classical spini^ 5 * 6 . While this approximation has been 
criticized^ 7 * 8 * 9 - it has been shown to yield good results at 
low but finite temperatures^ As the ferromagnetic cou- 
pling Jh of the itinerant electrons to the corespins is far 
larger than the hopping strength or the antiferromagnetic 
superexchange J', energy scales can be separated and the 
model can thereby be considerably simplified. In the cus- 
tomary Jh — ► oo approach and in the second order per- 



turbative treatment of the virtual excitations^ - for large 
but finite Jh, it is assumed that the e g electrons are par- 
allel to the local corespin. Double occupancies are thus 
suppressed and Coulomb repulsion is usually neglected. 

The FM Rondo lattice model with clas- 
sical corespins has been extensively investi- 
gated, 5 . 6 . 10 i n . 12 i 13 . 14 i 15 i 16 i 17 i 18 i 19 The influence of 
classical phonons^SiSi*^ and disorder* 2 **^ has been ad- 
dressed. Two-orbital models have been treated, mainly 
with classical phonons but without Coulomb repulsion. 
Aliaga et al. have examined the two-dimensional one- 
orbital Rondo model for Jh — > oo with an algorithm 
similar to oursi 9 " and reported several phases (stripes, 
island phases for commensurate fillings, and the so called 
"Flux Phase") as well as phase separation (PS), which 
has also been related by other authors ^i^i^&Sk Stripes 
were reported in a Rondo model applied to cuprates 
with different parameter values^SiSL 2 ^ 

In previous papers ji ** 7 * 1 ^ we have studied the one- 
dimensional model by Monte Carlo techniques similar to 
those used before. Careful analysis of the data showed, 
however, that instead of phase separation there are inde- 
pendent small FM polarons. We have recently extended 
these studies to the two-dimensional caseS 9 - for doping 
< x < 12% and at AFM corespin-coupling J' = 0.02. A 
thorough examination of a variety of observables showed 
that, like in ID, there is no phase separation but holes 
enter the system by forming independent FM polarons 
with a single hole inside. 

In the present paper we provide a phase diagram of 
the two-dimensional model as a function of corespin cou- 
pling J' and hole doping up to x — 0.6. We thoroughly 
examine the transition from the polaronic to the homoge- 
neous ferromagnetic phase. At moderate to large J' the 
polarons persist, whereas at very small J', phase separa- 
tion occurs. In the homogeneous phase, strong FM order 
is only present at small to moderate J'. 

At least for small doping, polarons are favored over 
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phase separation by entropy, and their dependence on 
temperature is therefore of interest. Previous studies 
have mostly employed an inverse temperature of (3 = 50 
(corresponding to an experimentally relevant tempera- 
ture of 50-100K, depending on hopping strength). We 
show that further lowering the temperature actually 
strengthens the polarons, even though the effect of en- 
tropy is reduced, i.e., they are also energetically favored. 
Only at small J', low temperatures enhance a tendency 
to phase separation and toward diagonal chains of holes. 

We also investigate the 3D model where we likewise 
find independent polarons at small doping, and in ad- 
dition a 'polaron lattice' at higher doping. No phase 
separation is observed at any doping. 

The outline of this paper is as follows. In Sec. |n] wc 
introduce the Hamiltonian and the numerical method. 
In Sec. IIIII we establish FM model polarons in two and 
three dimensions; we give an analytic estimate of their 
energy compared to phase separation in Sec. IIII Al and 
find polarons to dominate except for very small J' . Sec- 
tion If VI contains unbiased MC results for the two dimen- 
sional model. We present a phase diagram for (3 = 50 in 
Sec. UVDl 

We study temperature effects in Sec. II V El and a phase 
diagram for lower temperature (3 = 80 (30X — 60K) is 
presented in Sec. IIV E 41 

Finally, we present unbiased MC data on the three 
dimensional model at < x < 0.4 in Sec. W\ 

II. MODEL HAMILTONIAN AND METHOD 

In this paper, we treat the ferromagnetic Kondo lat- 
tice model in two and three dimensions with one orbital 
and classical corespins. We treat configurations with the 
electron spin antiparallel to the local corespin in second 
order perturbation theory, as proposed in Ref. IToL which 
is a systematic improvement over the Jh — > oo approach. 
The resulting effective spinless fermion (ESF) Hamilto- 
nian is 

<i,3> i,j <id> 

(!) 

The spinless fermion operator Cj (c\ ) destroys (creates) a 
local spin-up electron (i.e. an electron with spin parallel 
to the corespin at the same site) at site i. Down electrons 
have been integrated out, and the spin index has there- 
fore been omitted. Since the local corespins may point in 
an arbitrary direction, the ESF model Q still contains 
contributions from both spin orientations with respect to 
a global spin-quantization axis. Jh is the large ferromag- 
netic coupling of the e g electron spin to the corespin and 
is set to Jh = 6 throughout this work. The small AFM 
superexchange parameter J' takes values between J' = 
and J' — 0.05 for manganites. 

The hopping strength tjj in the kinetic energy (first 
term of Eq^l depends on the relative angle between 



the corespins at the two neighboring sites: 



= t (c iCj + SiSjjfa-M) = t Q cob(0«/2) , ( 2 ) 

where ct = cos(#j/2) and Sj = sin(??i/2). The fac- 
tor depends on the polar coordinates § and <f> of the 
corespins and contains a complex phase factor, which 
can lead among o thers to the so called "Flux Phase", 
see Refs. Il9l3fl3ll and Sec. WU We use the parameter 
to = 1 as unit of energy throughout this paper. 

The second term in Eq. ^ contains the second order 
treatment of the virtual excitations, where an up-electron 
at site i can become a down-electron at site j and imme- 
diately hop back to become an up-electron at site i again. 
Its strength also depends on the corespins but is always 
real: 



= toQdSje-^ + Si^e-^f) = i sin(^/2) 2 . 

(3) 

We treat this Hamiltonian by unbiased Monte Carlo sim- 
ulations for the grand canonical ensemble using the stan- 
dard algorithm- and for the canonical ensemble with the 
algorithm proposed in Ref. l29i The lattice size is 12 x 14 
sites in 2D and 6 x 6 x 4 in 3D where not otherwise indi- 
cated. We mainly choose non-quadratic lattices in order 
to minimize finite size effects due to closed shells. For our 
MC results, 50 to 200 lattice sweeps have been skipped 
between measurements and autocorrelations have been 
analyzed. Wherever they are found to be very long (no- 
tably around the phase boundary between the Flux phase 
and the FM/PM phase and for the vertical stripe phase), 
parallel tempering22i22i2i has been employed. 

III. FM MODEL-POLARONS IN 2D AND 3D 

In Ref. 12^. we showed that the results of unbiased 
Monte Carlo simulations for J' = 0.02, (3 = 50 and 
,/h = 6 for doping levels up to x ss 12% (i.e. near the 
completely filled lower Kondo band) correspond to inde- 
pendent ferromagnetic polarons and not to phase separa- 
tion. To this end, we constructed a model of independent 
polarons, evaluated several observables (spectral density, 
DOS, corespin correlation and dressed corespin correla- 
tion Eq. ©) for this model, compared the results to the 
MC data and found almost perfect agreement. 

In the undoped case x = 0, the corespins are aligned 
antiferromagnetically because of virtual excitations into 
the upper Kondo band and because of the AFM superex- 
change J'. The effective strength of the antiferromag- 
netic interaction is given by J c s = J' + jj-^ f° r n ~ I- 10 
One FM polaron in 2D consists of just one spin flipped 
from the AFM background and thus forming a fiv e-site 
FM region where a hole is trapped, see Fig. |l(a)| The 
eigenstates of the hole inside this small FM cluster lead 
to signals at w = ±2 and at uj = in the spectral density 
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FIG. 1: Idealized two-dimensional FM polaron of Lf = 5 
lattice sites, embedded in an AFM background (shown in 
Ref. I29I) (a) Spin and hole-density configuration for the 
groundstate. Empty (filled) squares represent spin down (up). 
Height represents hole density, (b) Contribution of the po- 
laron to the one-particle spectral function. For visibility, 
the 5-peaks in the spectral density have been broadened to 
a width of 0.2. 




FIG. 2: Spin and hole-density configuration of an idealized 
three-dimensional FM polaron of Lf = 7 lattice sites, embed- 
ded in an AFM background. The three layers represent three 
consecutive plains in the 3D lattice. Empty (filled) squares 
represent spin down (up). Height represents hole density in 
the ground state. 



and the DOS, see Fig. |l(b) For a large number of po- 
larons, i.e., of nipped spins, the FM polarons occasionally 
overlap and thereby form larger FM areas; such overlaps 
are neither suppressed nor encouraged in the independent 
polaron model. 

In 3D, the situation is very similar: When one spin 
is flipped in the AFM lattice, a FM domain of seven 
sites results. This is schematically depicted in Fig. [5] Its 
eigenenergies are given by to — (fourfold degenerate) 
and u) = ±v6- When a hole is inserted into the 3D 
lattice, it can gain the energy — V6, but the energy 12 J e s 
for the breaking of six AFM bonds has to be paid. The 
polaron energy is therefore given by 



e po i = — v6 + 12 J< 



eff 



(4) 
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FIG. 3: (Color online) Critical chemical potential for phase 
separa tion ( solid line) and FM polarons of the shape depicted 
in Fig. |l(a)|for (a): 2 dimension and (b) 3 dimensions. 



This energy determines the critical chemical potential 
fx* = — e po i, at which holes that form polarons enter the 
completely filled lower Kondo band of a 3D system. In- 
terestingly, such seven-site stars were also reported to 
result from combined lattice and spin effects for the elec- 
tron doped system, i.e., for few electrons^ 



A. Energy comparison of Phase Separation and 
FM Polarons 

We examine the competition of phase separation and 
FM polarons by comparing their energies for the relevant 
range of the effective antiferromagnetic superexchangc 
Jeff = J' + 277/7 which can of course only be done ap- 
proximately. We obtain the critical chemical potential 
/ips for phase separation by setting the energy of the 
completely filled AFM band equal to that of a partially 



filled FM band. 

/MPS rf-ps 
e n(e) de~[ips / n(e) de , (5) 
-z J — z 

ekin 

where z is the number of nearest neighbors for each site 
(4 in 2D, 6 in 3D) and n(e) is the one-particle den- 
sity of states of a tight binding band with e{k x ,k v ) = 
—2cos(k x ) — 2cos(k y ) or e(k x ,k y ,k z ) = —2cos(k x ) — 
2 cos(ky) — 2 cos(fc z ), respectively. The kinetic energy per 
site ekin (first term on the right-hand side) is zero for the 
completely filled band. 

The critical chemical potential for FM polarons is given 
by Mpoi = 2 — 8 J c ff in 2D (see Ref. l29^ and ^ po i = — 
12J e ff in 3D (see EqEJ. It is compared to fip$ in Fig.El 
for two and three dimensions. In contrast to the energy 
minimization performed for ID systems (given in Ref. llq) 



4 



we did not optimize the shape of the polarons, but we 
always used polaron s con sisting of one single flipped spin 
as depicted in Figs. 1(a) and [5] This puts FM polarons 
slightly at a disadvantage. 

For \x larger than the critical values, the band is com- 
pletely filled and AFM. When the chemical potential is 
lowered and holes are introduced into the system, this 
leads to PS if /ips > M po i and to polarons if fi po i > /ips- 

For 2D, one sees that PS is favored for small J c g < 0.08 
in accordance with the MC data, see Sec. IIVBI For the 
parameter Jh = 6 used in our calculations, J c ff sw 0.083 
corresponds to J' = 0, where we do indeed observe phase 
separation. For larger values of J c g, i.e., smaller Jh < 6 
or larger J' > 0, polarons are not only favored by entropy, 
but also by energy. 

In three dimensions, the preference for polarons over 
PS sets in for even smaller values for J e s- 



IV. MC RESULTS FOR 2D 

We first investigate the influence of AFM superex- 
change J' and doping x at f3 = 50, eventually leading 
to the phase diagram in Fig. Temperature effects are 
discussed in Sec. IIV El with a phase diagram for /3 = 80 
in Fig. EH 




(b) 





FIG. 4: MC snapshot of the hole density for (a) 20 holes 
(x » 12%, shown in Ref. 2Sj), (b) 24 (x » 14%), and (c) 31 
(x m 18%) holes in a 14 x 12 lattice at 13 = 50, J' = 0.02, J H = 
6. Height represents hole density, grayshades are for better 
visibility. 




A. J' = 0.02: From independent polarons to phase 
separation 

In Refi ll we found a very close correspondence of the 
MC data with the results from the polaron model for all 
observables up to doping levels of x s=s 12%, where the 
polarons already cover approximately 60% of the lattice. 
The polarons appear to be independent, with overlapping 
polarons occurring at a similar rate as for the indepen- 
dent polaron model. 

In this section, we investigate higher doping levels 
x > 12% at /3 = 50, where the results begin to deviate 
from the independent polaron results; the homogeneous 
FM phase sets in at x w 21%. In between, the polarons 
attract each other which leads to phase separation (PS). 
The transition from polarons to PS is not well defined. 
The polarons first coexist with large clusters, which fi- 
nally dominate. This development is illustrated in the 
MC snapshots depicted in Fig. 0] It shows snapshots 
for 20 holes (x » 12%), where 19 polarons can be seen 
and only one hole is delocalized, for 24 holes (x ~ 14%), 
where polarons coexist with a larger and more homoge- 
neous area, and finally for 31 holes (x w 18.5%). This 
last doping is only approximately 4 holes away form the 
homogeneous FM phase at x > 21%, and even there, 
some polarons persist. 

In order to check whether the addition of holes leads 
primarily to more small polarons or to a growth of the 
existing ones, we compute a dressed corespin correlation 



FIG. 5: (Color online) Dressed corespin correlation Eq. © for 
24 holes (x, dashed), 31 holes (o, dotted) and » 36 (A, dash- 
dotted) holes. The inset shows the usual corespin correlation 
in real space S(r) = j- Remaining parameters as 

Fig. [I| 

function 

i 

The hole density at site i, denoted by ni, is related to 
the electron density via nh = 1 — = 1 — (etc-) 5, where 
( n i)s is the expectation value of the density given the 
corespin configuration S. The sum over i is taken over 
all lattice sites and the observable is averaged over all 
corespin configurations S occurring in the MC run. The 
dressed correlation measures the ferromagnetic regions 
around holes. 

Figure [S] shows Sh(r) for 24, 31 and w 36 holes. The 
latter doping marks the point, where the more homo- 
geneous phase sets in and where the compressibility be- 
comes much smaller than in the polaronic and phase sep- 
arated regimes, so that this filling can be obtained in the 
grand canonical ensemble. One sees that the ferromag- 
netic regions around the holes grow slightly with doping, 
e. g., the correlation at r 2 = 5 goes from AFM to FM. 
This is an indication for phase separation. The usual 
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corespin correlation in real space is shown in the inset. 
It shows the evolution of short-range ferromagnetic cor- 
relations. However, no long-range FM correlations are 
present even for sw 36 holes (x s» 21%). This shows that 
the orientation of the small FM clusters is largely inde- 
pendent, see also the spin structure factor in Sec. IIV CI 
We therefore do not have coexistence of macroscopic FM 
and AFM phases at any point for J' = 0.02. 



B. Dependence of Polarons and Phase Separation 
on J' 



Although J' is usually taken to be very small com- 
pared to the other parameters in the Hamiltonian, it has 
a considerable influence on the polaronic regime, because 
it stabilizes the AFM background. 



(a) 



(b) 





FIG. 6: Monte Carlo snapshots for J' = and /3 = 50. (a) 
6 holes, polarons; (b) 20 holes, phase separation. The MC 
simulations were done for Jh = 6, /3 — 50 on a 12 x 14 lat- 
tice. Height represents hole density, grayshades are for better 
visibility. 

For J' = and f3 = 50, MC snapshots shown in Fig. 
reveal polarons with a tendency to form clusters for a few 
holes and phase separation for larger doping. 




FIG. 7: (Color online) Dressed corespin correlation Eq. © 
from unbiased MC Data for J' = (symbols) and J' = 0.02 
(continuous lines) for 1 (o, dotted), 6 (x, solid) and 20 (V, 
dashed) holes. The inset shows the corespin correlation. Pa- 
rameters as in Fig. 

The dressed spin correlation Eq. |J5J| for 1, 6 and 20 
holes is depicted in Fig. □ for J' = and J' = 0.02. In 





— 3 holes, x = 2% 


(a) A 


---6 holes, x « 3.5°/ 


12 holes, x = 7% 




■—■20 holes, x=12% 








FIG. 8: One-particle DOS for (a) J' = 0.02 (shown in Ref. 
I29I) and for (b) J' = 0, other parameters as in Fig. OH 



this doping range, the ferromagnetic regions around the 
holes do not grow with doping for J' — 0.02 (continuous 
lines), instead, the nearest AFM correlation at r 2 = 5 
remains AFM, which indicates polarons. In contrast, the 
FM regions do grow for J' = (symbols), this growth 
suggests larger FM regions, i.e. phase separation. Also, 
the antifcrromagnetic correlations shown in the inset de- 
crease faster for J' = than for J' = 0.02. 

Fig. shows the DOS for J' = 0.02 and for J' = 
at various fillings. For both values of J' and all de- 
picted doping ranges, one sees the central band which 
stems from the movement of the electrons in the imper- 
fect AFM background. This band is completely filled. 
For small doping, one sees the polaronic states atwsj ±2, 
see also Fig. ^ For J' = 0.02, the weight of these peaks 
increases with doping, but their shape changes only for 
very large doping, in accordance with independent po- 
laron results^ The pseudogap remains clearly visible. 
For J' = 0, one also sees the polaronic states and the 
pseudogap for small doping, but the pseudogap is filled 
at larger doping. The states filling the pseudogap stem 
from larger FM clusters, and these are larger and more 
frequent than can be accounted for by occasionally over- 
lapping independent polarons. The actual experimental 
observation of the pseudogap^iiSLSiSi rather favors the 
polaron scenario as opposed to PS. 

Figure [5] shows the spectral density for J' = 0, J' = 
0.02, and J' = 0.05 for 20 or 21 holes, i.e. x m 12%. 
The larger FM areas for J' = [top panel (a)] lead to a 
second band in addition to the one stemming from the 
AFM parts of the lattice. It has a larger bandwidth, be- 
cause the corespin correlations are ferromagnetic and the 
hopping strength is therefore larger. The central panel 
(b) shows the results for the independent polarons at 
J' = 0.02 (also given in Ref. 12^) . Around the polaronic 
states at uj = ±2, one sees the signals from overlapping 
polarons, which are also there for the independent po- 
laron models The signals from the small FM clusters 
do not form a continuous band, but are instead separated 
from the AFM band by the pseudogap. 

A stronger antiferromagnetic superexchange J' = 0.05 
suppresses configurations with overlapping polarons, as 
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(a) 3 



J'=0 




(b) 
J'=0.02 




(c) 
J'=0.05 




FIG. 9: Spectral density for x « 12% and (a) J' = (20 
holes), (b) J' = 0.02 (20 holes, also shown in Ref. H|, and 
(c) J" = 0.05 (21 holes) on a 12 x 14 lattice, f3 = 50, J H = 6. 



can be seen in the lowest panel, Fig. |9(c)| In addition 
to the antiferromagnet, the data for J' — 0.05 show al- 
most only one-polaron states and the pseudogap is more 
pronounced. This shows that not even overlapping po- 
larons exist for J' = 0.05, let alone phase separation. 
The reason for this behavior is that larger J' favors an- 
tiferromagnctically stacked polarons over larger clusters. 
This can be seen in another dressed corespin correlation, 
which takes into account hole densities at both lattice 
sites: 



i+r i+f 



(7) 



Note that even without the corespins, this would not 
give the usual density correlation, because it measures 
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ally, this observable is of course averaged over all spin 



FIG. 10: (Color online) Hole-spin-hole-spin correlation as de- 
fined in Eq|7| for 20 holes and various values of J' : Other 
parameters as in Fig. El 



configurations S obtained in the Monte Carlo run. This 
correlation is depicted in Fig. ^| for 20 holes (x « 12%) 
and for J' = 0, J' = 0.02 and J' = 0.05. One sees 
a positive correlation indicating PS for J' — 0, almost 
only short-range correlations over two sites (one polaron) 
for J' = 0.02 and a negative correlation at r 2 = 5 for 
J' = 0.05, which indicates antiferromagnetically stacked 
neighboring polarons. 



C. FM, PM and Flux Phase 

The exchange interaction also plays an important role 
at electron densities below the polaronic and phase sepa- 
rated regimes, i. e., at larger doping levels x > 0.2, which 
corresponds to chemical potentials jj, < fi* = — e po i. 

For small to medium J' < 0.03, this region is ferro- 
magnetic because of the dominant double-exchange. The 
corespin structure factor S(k) is depicted in fig. 1111 for 
J' = 0.02 and variable filling. The curve for x « 12%, 
which is in the polaronic regime, shows a clear AFM sig- 
nal. There is a barely visible FM signal, but it is very 
small in spite of the rather large doping, because the po- 
larons are independent and not mutually aligned. Just 
below fj,*, corresponding to x w 21%, the AFM peak 
has almost vanished and there is a slightly larger FM 
peak at (0,0). This means that the ferromagnetism is 
not very strong and that even some remnant of the AFM 
remains, see also the corespin correlation in real space 
in Fig. The lattice is however homogeneous and the 
spectral density (not shown) consists of one single tight 
binding band with reduced hopping strength t » 0.75 
and with broadened signals. For larger doping x « 29% 
and x w 32%, the spin structure factor S(k) shows that 
the ferromagnetism grows upon increasing the hole den- 
sity up towards half-filling of the lower Kondo band. The 
bandwidth then increases and the signals in the spectral 
density become sharper, as fluctuations of the corespins 
are suppressed by the kinetic energy. 
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20 holes = 12 % 
* 36 holes = 21 % 
« 48 holes = 29 % 
= 54 holes = 32 % 



(0,0) 



(0,7t) (JI.JI) 

wave number 



(0,0) 



FIG. 11: Spin-spin correlation in momentum space for various 
hole doping, J' = 0.02, J H = 6 and j3 = 50. 



With increasing J' , however, the ferromagnetism grad- 
ually becomes weaker and for J' = 0.05, S(k) does not 
show any signal at (0, 0) for any doping x < 0.5. Instead, 
the so called "Flux phase" appears around half filling. 
The band structure then has a pseudogap at half-filling 
and differs markedly from the tight-binding-DOS of a fer- 
romagnetic or paramagnetic lattice. The spin structure 
factor shows signals at (0,7r) and (it, 0)i 19 i 30 i 31 i 40 



D. Phase Diagram for — 50 

In order to determine the phase diagram, shown in 
Fig. 1121 criteria to distinguish the phases have to be 
specified. For the crossover from the polaronic to the 
phase separated regime, we choose the filling, at which 
the ferromagnetic regions around the holes - detected by 
dressed corespin correlation Eq. © - begin to extend, 
i. e., when the nearest antiferromagnetic correlations at 
r 2 = 5 become ferromagnetic. It must be emphasized 
that the transition is a gradual one: larger ferromagnetic 
domains occasionally occur at smaller doping and above 
all polarons persist to larger doping (Sec. II V A[) . This 
criterion is therefore somewhat arbitrary. 

To determine the phase boundary between the pola- 
ronic/phase separated region and ferromagnetism on one 
hand, and between the polaronic regime and paramag- 
netism on the other, we choose the electron density just 
below the discontinuity at the critical /x*, although the 
antiferromagnetic peak may still be visible at this point; 
the distinction is therefore rather based on the compress- 
ibility. From there on, a system is labeled "ferromag- 
netic" or "paramagnetic". 

In the paramagnetic regime, we find hardly any mag- 
netic structure for intermediate doping, until the signals 
from the Flux phase begin to appear. The crossover 
from FM to PM is continuous and depends on doping, 
larger x leading to more ferromagnetism. As a criterion 
for the transition from FM to Flux phase, one can take 
the maximal change of the corespin moment with filling 
d\Stot\/dN , because |5t Q t| = for the ideal Flux phase. 
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FIG. 12: Phase diagram for J H = 6, = 50, < J' < 0.05 
and < x < 0.6 (i.e. filled to 40% filled lower Kondo band) 
on a 14 x 12-lattice. "Pol.": polaronic regime, "Pol. /PS": 
mixture of both polarons and larger ferromagnetic clusters, 
"AFM": antiferromagnet, "FM": ferromagnet, "PM": regime 
without magnetic structure, "Flux" : Flux phase. 



Obviously, this becomes less exact as the ferromagnetism 
is weakened. A further indicator is the change in the 
spectral density, which for the Flux phase differs very 
much from a tight-binding one. The Flux phase is sub- 
ject to strong finite size effects, which are caused by the 
fact that only a very small part of the Brillouin zone con- 
tributes to the few states around the pseudogap around 
(7r/2,7r/2). For this reason, lattices of different sizes, 
above all square systems because of closed shell effects, 
may exhibit the transition from ferro- /paramagnetism to 
the Flux phase at slightly different filling. For hole dop- 
ings above the Flux phase, all systems become ferromag- 
netic. 

The resulting phase diagram with the phases discussed 
above is shown in Fig. 1121 



E. Effects of a lower temperature 

1. J 1 = 0.02: Polarons stabilized further 

For J' = 0.02, where the system can be well described 
by independent polarons at (3 = 50 and x < 12%, cool- 
ing to P = 80 suppresses overlapping polarons, larger 
ferromagnetic regions and phase separation. The sup- 
pression of overlapping polarons can, e.g., be seen by 
comparing the spectral density at doping x — 12% (20 
holes) for J' = 0.0 2 and (3 = 80 (Fig. [T3J| to the one for 
(3 = 50 (Fig. |9(bJ |. Besides the narrowed antiferromag- 
netic band, almost only one polaron signals are seen at 
the lower temperature, which means that the polarons 
are not independent but avoid overlapping. 

The suppression of phase separation by lower temper- 
ature is clearly visible in the MC snapshots for 28 holes 
(x ~ 17%) depicted in Fig. [21 While only a few polarons 
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FIG. 13: Spectral density (same as fig. |9(b)| but for lower 
temperature f3 = 80.) 




(a) (b) (c) 



FIG. 14: MC snapshot of the hole density for 28 holes (x « 
17%) at temperatures (a) /3 = 50, (b) (3 = 80, (c) and /3 = 100. 
Other parameters as in Fig. [I] 

exist for (3 = 50, they dominate at (3 = 80, where only a 
few holes are delocalized, and even more so at (3 = 100, 
where two polarons hardly ever share a site. 

As for (3 — 50 and J' = 0.05, the reason for the 
polarons' dominance over PS is their antifcrromagnetic 
stacking. Since each polaron takes 5 lattice sites, they 
can obviously no longer be stacked in this way for dopings 
larger than 20%. As can be seen in the phase diagram in 
Fig. this doping is the point, where the larger clusters 
begin to dominate. In this case, the polaronic phase is 
more ordered than a phase separated scenario because of 
the stacking, and PS therefore occurs at higher temper- 
ature because of its higher entropy. 

The stacking can be seen in the Hole-Spin-Holc-Spin 
correlation depicted in Fig. 1151 for various doping levels. 
While one sees only short range correlations for 6 holes 
(x = 3.5%), the antiferromagnetic stacking can be clearly 
observed for 28 holes (x = 17%). This order disappears 
again for 34 holes (x — 20%), when the polarons begin 
to merge. At this point, some tendencies to diagonal 
chains (see below) may exist at very low temperatures 
03 = 100). 

In the ferromagnetic part of the phase diagram, a spin 
canted phase becomes actually more stable than ferro- 
magnetism for some electron fillings and for periodic 
boundary conditions. Its spin structure factor has a peak 
for the smallest reciprocal lattice-vector of the system 
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FIG. 15: (Color online) Hole-Spin-Hole-Spin correlation Eq0 
for J' = 0.02, (3 = 80 and various fillings. Other parameters 
as in Fig. 
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FIG. 16: Corespin configuration with t hree flipped spins, 
filled with two holes. Details as Fig. |l(a)| 

(ir/7 for a 14 x 12 lattice). Analytic energy comparison of 
the FM and this canted phase shows that the range of fill- 
ings, where the canted phase has lower energy, becomes 
smaller for larger lattices. Moreover, it is not present 
for open boundary conditions and thus probably a finite 
size effect. It shows, however, that ferromagnetism is not 
very dominant for doping levels near the phase separated 
regime. 



2. J' — 0: Tendency to diagonal stripes 

In Ref . |2(J we showed that larger FM clusters with two 
flipped spins and two holes have (slightly) higher energy 
than two separated polarons. Fig. ^j] depicts another 
possibility with three flipped spins but only two holes. 
The total kinetic energy of the two holes is then 0.61 to 
lower than that for two independent polarons. However, 
3 spins have to be flipped instead of two, which costs the 
energy 8 x J e g. The gain and the loss balance each other 
for J c ff w 0.077 in this simple model calculation. In the 
MC simulations, we observe tendencies to such configura- 
tions for J' = (J c ff w 0.083), which extend to diagonal 
chains at higher doping, see the snapshots in Fig. For 
larger hole doping, isotropic clusters coexist with these 



(a) 



(b) 



(a) 



(b) 



FIG. 17: Monte Carlo snapshots for J' = 0.0 and (3 = 80 
and (a) 6 holes (x « 4%), (b) 12 holes (ac ~ 7%). Height 
represents hole density, grayshades are for better visibility. 





FIG. 18: (Color online) Hole correlation Eq. @ for J' = 
0.0 and 12 holes (x = 7.1%): The left panel (a) shows the 
correlations in (1, 0)-direction; the right panel (b) shows the 
same for the diagonal (1, Indirection. Dashed line: (5 = 50, 
solid: = 80. 



diagonal chains and finally dominate. In order to mea- 
sure these diagonal chains, we use the hole-correlation 
function 



-n h ) (n^-ffc), (8) 



where n? is the expectation value of the hole density at 

site i for a given corespin configuration and is the 
overall hole density. As explained for Eq. this is not 
the usual density correlation function. The results for 12 
holes (x = 7.1%) are depicted in Fig. 1181 and show that 
diagonal correlations are slightly enhanced at {3 = 80. 



3. J' — 0.05: Vertical stripes 

For J' = 0.05 and doping x = 0.25, low tempera- 
tures favor a stripe-phase with a periodicity of 4 sites 
perpendicular to the stripes. The spin configuration is 
schematically depicted in Fig. 
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The corespins are 
drawn within the xy-plane for visibility, but such config- 
urations are not preferred over others with the corespins 



FIG. 19: (Color online) Stripe Phase for J' = 0.05 and 
x = 0.25. Left: Schematic representation of the corespins. 
The circles represent the holes density, but only qualitatively. 
Right: one-particle density of states of the two dispersive 
bands for Jh — * oo and the thermodynamic limit (contin- 
uous line) and for finite Jh on a 12 x 14 lattice with random 
fluctuations added to the corespins. 



perpendicular to the plane of the 2D-lattice, because the 
Hamiltonian only depends on the relative orientation of 
neighboring spins, and a global rotation therefore has no 
effect. 

For the schematic corespin configuration with perfectly 
aligned spins and with the approximation Jh — > oo, the 
band structure is gapped and consists of several disper- 
sionless bands at u> — and two bands with e(k x ) = 
±y/4cos 2 (k x ) 4- 2, when the stripes run in x-direction. 
As no hopping is possible between the antiparallel stripes, 
the band does not depend on the momentum perpendic- 
ular to the stripes. The resulting one-particle density of 
states for Jh — > oo and the thermodynamic limit can be 



seen in Fig. 19(b) (continuous line) and it is very similar 
to the DOS of one-dimensional tight-bin ding ba nds. For 
comparison to the MC data (inset of Fig. [20(b) I, we also 
calculated the spectral density for finite Jh on a 12 x 14 
lattice for the model corespin configuration with slight 
additional fluctuations (dashed line in Fig. 19(b) I. The 
fluctuations were chosen so as to yield approximately the 
same width for the central band as the MC data. On a 
24 x 12 lattice, the stripes can only develop in one direc- 
tion, because 14 is not divisible by 4. One therefore can 
see the dependence of these bands on only one direction 
of the momentum quite clearly in t he ave raged spectral 
density depicted in Fig. 20(b) Fig. 20(a) shows an MC 
snapshot of this phase, where one sees stripes in the hole 
density. 

Figure shows the Hole-Spin-Hole-Spin correlation 
function Eq[7|for x — 0.25: One clearly sees the regular 
array of holes and the stripes in the corespins as depicted 



in Fig. 19(a) 



The compressibility around this filling is very large, 
i.e., there is a jump in the electron density versus the 
chemical potential from x ~ 0.29 toisj 0.23 and the fill- 
ing x — 0.25 can not be stabilized in the grand canonical 
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(a) (b) 

FIG. 20: MC snapshot and spectral density for J' = 0.05, fi = 
80, and x — 0.25. Left: MC snapshot: The heights repre- 
sent hole densities, grayshades indicate FM (white) and AFM 
(black) correlations to the nearest neighbors. Right: Spectral 
density and one-particle DOS. 




1 to stripes 



FIG. 21: Hole-Spin-Hole-Spin correlation Eqd for J' = 
0.05,/? = 80 and x = 0.25 (42 holes). Grayshades for bet- 
ter visibility. 



ensemble. The discontinuity is not present for the higher 
temperature (3 = 50, where there is no stripe phase. At 
fillings slightly away from x — 0.25, stripes develop in 
parts of the lattice and are mixed with polarons or larger 
clusters. 

The occurrence of this stripe phase at x — 0.25 and 
J' = 0.05 is in contrast to the results reported by Aliaga 
et ali^, who found an island phase of antiferromagnet- 
ically stacked FM clusters of size 2 by 2 for this filling. 
The electron density is homogeneous in this phase, and 
every spin has two FM and two AFM neighbors. Their 
calculations were done for the Jh — > oo model, but with 
a larger J', so as to give approximately the same effec- 
tive parameter J e g. This island phase can however not be 
stable on a 12 x 14 lattice. In order to check the validity 



FIG. 22: Phase diagram for /3 = 80 and (3 = 50. Similar as 
fig. 1121 but for a smaller range of doping. Solid lines: results 
for P = 80, dashed lines: (3 = 50. 

of our results, we therefore repeated the simulations on a 
12 x 12 system and still found stripes. Numeric compar- 
ison of the ground-state energies for perfect stripe and 
perfect island phase on a 12 x 12 lattice gave a lower 
energy for stripes for Jh = 6 and J' = 0.05 (correspond- 
ing to Jeff w 0.133 for x = and to J c s 0.1125 for 
x = 0.25), but favored the island phase for Jh = oo 
and J' = 0.133 or J' = 0.1125. The reason lies in the 
density dependence of the AFM coupling coming from 
the virtual excitations (second term in Eq. , because 
the sites with 3 AFM and only one FM neighbor have 
high electron density and therefore maximize the gain 
from this interaction, whereas the holes weaken the in- 
teraction on the sites with predominately FM neighbors. 
This therefore presents a feature which is missing from 
the first order Jh — > oo approximation but is included in 
the second order approach to the FM Kondo model. 

4- Phase diagram at (3 = 80 

Our findings are summarized in the phase diagram for 
/3 = 80, < x < 0.36 and < J' < 0.05 depicted in 
Fig. [52 It shows that lower temperatures lead to larger 
clusters and PS for small J' < 0.01, while they favor 
individual polarons for more realistic values J' > 0.01. 
For the lower temperature, the homogeneous ferro-/para- 
magnetic phase only sets in at higher doping. 

V. MC RESULTS FOR 3D 

In three dimensions, we use a lattice of 6 x 6 x 4 sites 
with periodic boundary conditions. For selected doping 
levels, we have repeated the simulations on lattices of dif- 
ferent size in order to control finite size effects. As in one 
and two dimensions, the corespins are antiferromagneti- 
cally aligned for the filled lower Kondo band and the sys- 
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FIG. 23: MC snapshots for (a) 3 holes (x « 2%), (b) 14 holes 
(a; ~ 10%), (c) 32 holes (x « 22%) and (d) » 38 holes (a; » 
26%). The four plots stacked one above the other represent 
the 4 layers of the 6x6x4 lattice, height represents hole 
density, grayshades are for better visibility. Jh — 6, /3 = 
50, J' = 0.02. 



tem exhibits a transition to ferromagnetism upon doping. 
Figure BBI shows MC snapshots for J' = 0.02 and = 50 
at various doping levels. For three holes (x ~ 2%), one 
clearly sees three polarons of 7 sites each. The polarons 
also remain well separated for 14 holes (x « 10%), al- 
though they take up w 68% of the lattice at this doping 
level. As every polaron contains 7 sites, no more than 
x w 14% holes can be doped into the lattice in this way. 
For the larger doping level x ss 22% (32 holes), one sees 
tendencies toward a superstructure of regularly arrayed 
holes, see Fig, ^f c). which we termed 'polaron lattice'. 
An idealized representation of the polaron lattice is de- 
picted in figure 1241 which is for a doping concentration 
of x = 0.25. The idealized version consists of two sublat- 
tices with opposing corespin and hole density is high on 
sites with four FM neighbors. Interestingly, at this 'op- 
timal' doping concentration for the polaron lattice, the 
actual MC simulations favor a more homogeneous and 
disordered structure, like the one depicted in Fig. EJd), 
which has been obtained for x as 27% (ss 39 holes). 

On very small systems (4x4x4 and 4x4x6), 
the polaron lattice competes with a three-dimensional 
form of the vertical stripes found in 2D for J' = 0.05 
and f3 = 80 (Sec. IIVE3|) . Due to their periodicity, the 
stripes can only develop on lattices, where at least two 
of the three dimensions are divisible by four, and they 
do therefore not occur on a 6 x 6 x 4-system. Numeric 
comparison of the energies of the polaron lattice and the 
stripe-structure for x = 0.25 on systems of different size 
revealed that the polaron lattice has lower energy for sys- 
tems which are larger than 4 in at least two directions, 
i.e., for lattices which are not of the shape 4 x 4 x L z . 




FIG. 24: Schematic representation of the corespins and den- 
sity for the 'polaron lattice'. Empty(filled) spins denote spin 
up (down), height represents hole density. 



In order to verify this result, we performed Monte Carlo 
simulations on an 8 x 6 x 4-system and indeed found a 
polaron lattice. 




M 



FIG. 25: Electron density n e \ — 1 — x vs. chemical potential 
/j, in 3D. Parameters as Fig. 1231 

The dependence of the electron density on the chemical 
potential is given in Fig. 123 The slope of the curve n e i 
vs. fi (i. e. the compressibility) is finite for the polaronic 
phase near n c \ = 1 in contrast to the results for inde- 
pendent polaronsji&Si because the polarons avoid over- 
lapping and it therefore becomes harder to fit more holes 
into the lattice. At x as 1/7 (n c \ as 0.86), when there is no 
more room for polarons, the situation changes. The slope 
of n c \ vs. fj, becomes steeper, and the polaron lattice par- 
tially develops until x < 23%, where the homogeneous 
disordered phase sets in. At x as 27%, the compress- 
ibility is suddenly much reduced. It is not clear, why 
the homogeneous phase has such a high compressibility 
for 0.23 < x < 0.27, and this may be a finite size ef- 
fect. The lattice is not phase separated in this regime, 
snapshots are similar to the one depicted in Fig. I23f d). 
with some sites (scattered throughout the lattice) having 
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much higher or lower density than others. After ferro- 
magnetism has become established at still higher doping 
values (x > 30%), the density becomes more evenly dis- 
tributed (not shown). 
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FIG. 26: (Color online) 3d Corespin Structure Factor for var- 
ious doping levels. Parameters as Fig. 1231 

Figure 1261 shows the momentum dependent corespin 
correlations for approximately the same dopings as the 
snapshots. For the almost filled band (3 holes, x w 2%), 
it shows antifcrromagnetism which decreases slightly for 
14 holes (x w 10%). For 32 holes (x w 22%), the partial 
polaron lattice, a second peak appears for k = (0,0, n). 
For the polaron lattice, according to figure l2*4l the major- 
ity spin for consecutive planes alternates, which explains 
the signal at (0,0, ir) in S(k). The corespin structure 
is isotropic, the same alternating structure occurs in x— 
and y— direction and analogous signals appear at (0, tt, 0) 
and (77, 0,0). For the more homogeneous phase at « 39 
holes there is no magnetic structure, and ferromagnetism 
only begins to set in for rj 45 holes (x rj 30%, not shown) 
and it is established for ss 58 holes (x rj 40%). 

Figure 12*71 shows the one-particle density of states for 
various doping levels. For 3 holes (x « 2%), one sees the 
broadened central band from the AFM background and 
the signals from the polaronic states at wi y6- For 14 
holes (x w 10%), these signals are stronger and some- 
what broader, because the polarons now cover a large 
part of the lattice and occasionally connect. For 32 
holes (x sa 22%), the polaronic states have a consid- 
erable width, but the pseudogap is well preserved. In 
will be shown later, that these bands are the dispersive 
bands of the 'polaron lattice'. The pseudogap has closed 
after the onset of the homogeneous phase at « 39 holes 
(x ri 27%), but some remnant of it is still visible. For still 
higher doping of rj 58 holes (x rj 40%), where the lattice 
is ferromagnetic, the band (not shown) is a tight-binding 
band with a slightly reduced hopping amplitude. 

In closing, let us discuss band structure and DOS for 
the ideal polaron lattice, depicted in figure |2"H The unit 
cell in this case contains eight sites, and for perfectly 
FM/AFM corespins, the one-particle DOS has a degen- 
erate non-dispersive band at u> = —2/Jh and two mo- 
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FIG. 27: (Color online) One-particle DOS for various doping 
levels. Parameters as in Fig. 1231 



mentum dependent bands at energies 



e(k x , ky, k z ) — — 



± ]JjT+ 4(cos 2 (fc a; ) + cos 2 (fc a ) + cos 2 (fc 2 )) . (9) 



The one-particle density of states resulting from this 
dispersion is depicted in Fig. [2H for a 6 x 6 x 4 lattice 
(dash-dotted line). The two dispersive bands consist of 
a number of peaks because of finite size effects. For an 
infinite lattice (solid line), the DOS has a pseudogap be- 
tween the dispersive bands and the central (5-peak. In ad- 
dition, we have performed a model calculation, in which 
the corespins of the ideal polaron lattice are randomly 
disturbed. The DOS of the resulting tight-binding model 
with transfer integrals given by Eq.|2]is compared to the 
result of the ideal polaron lattice in Fig- EHl The amount 
of fluctuations has been chosen so as to give approxi- 
mately the same wi dth for the central band as the MC 
data (similar to Fig. 19(b) I. 

While the model calculation (for x — 0.25) differs 
slightly from the MC data for 32 holes (a; RJ 22%) de- 
picted in Fig. [23 the width of the bands and their appar- 
ent double-peak structure is reproduced. The remaining 
differences are probably due to the fact that the MC data 
are for x = 0.22 and the polaron lattice is only partially 
developed, see the MC snapshot Fig. I^ST c). 

For J' = 0, diagonally arranged polarons similar to 
those observed in 2D (Sec. IIVE 2|) were found for x < 
0.25 on a 6 x 6 x 4 lattice. These stripes occurred on 
every second plane along the z~ direction (with L z =4). 
However, this phase was not observed on a 6 x 6 x 6 lattice 
and we therefore consider it a finite size effect. On this 
larger lattice, we observed the 'polaron lattice' instead 
and the results were generally very similar to those for 
J' = 0.02. 
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FIG. 28: (Color online) One particle DOS for the ideal po- 
laron lattice in the thermodynamic limit (solid) , for the ideal 
case on a 6 x 6 x 4 lattice (dash-dotted line) , and for 6x6x4 
with random fluctuations around the perfect spin arrange- 
ment (dashed line). 

VI. CONCLUSIONS 

We have used Monte Carlo simulations to examine the 
two- and three-dimensional ferromagnetic Kondo lattice 
model at small to medium hole doping and for parameters 
relevant to manganites. In three dimensions, we find that 
a small additional antiferromagnetic exchange J' = 0.02 
yields very similar results as J' = 0. No phase separation 
is observed at any doping level. Instead, there are again 
polarons consisting of a single flipped spin in an AFM 
background and containing a single hole, which form a 
regular 'polaron lattice' when they become very dense. 

In two dimensions, we have focused on the polaronic 
phase and its boundaries toward phase separation. At 
small hole doping, we find polarons for all values of the 
superexchange J'. For larger doping and small J', e.g., 



J' = instead of J' = 0.02, the polarons attract each 
other and tend to form larger clusters. Eventually, upon 
higher hole doping, phase separation into large ferromag- 
netic and antiferromagnetic regions sets in. A large an- 
tiferromagnetic superexchange (J' = 0.05), on the other 
hand, suppresses overlapping polarons and larger ferro- 
magnetic regions, because it stabilizes the AFM back- 
ground and favors antiferromagnetic stacking of the in- 
dividual polarons. In the doping range 0.2 < x < 0.5, J' 
suppresses the FM long-range order, so that the system 
is in reality rather paramagnetic than ferromagnetic for 
J' > 0.03. 

We find that lowering the temperature generally moves 
the homogeneous ferro-/para-magnetic phase to higher 
doping. Phase diagrams for = 50 and = 80 are given 
in Fig. ^] and Fig. |2U A lower temperature stabilizes 
the polarons and disfavors larger clusters for intermedi- 
ate and large J' > 0.01; for J' = 0.02, phase separation 
only occurs, when no more polarons can be fit into the 
lattice. PS is favored by lower temperatures only for un- 
physically small J' < 0.01. This is corroborated by a 
comparison of ground-state energies for idealized phase 
separation and polaron scenarios, which gives phase sep- 
aration for small J c g < 0.08 (corresponding to J' w 
for Jh = 6), and polarons for larger J c ff > 0.08. For 
J' w 0, lower temperatures enhance ferromagnetic di- 
agonal chains with delocalized holes, while they lead to 
vertical stripes for large J' = 0.05 at doping x = 0.25. 
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